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On the stability of the PWP method 
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Abstract 

The PWP method was introduced by Diaz in 2009 as a technique for measuring 
indirect influences in complex networks. It depends on a matrix D, provided by 
the user, called the matrix of direct influences, and on a positive real parameter A 
which is part of the method itself. We study changes in the method’s predictions 
as D and A vary. 


1 Introduction 

One of the main problems in network theory is to dehne a ranking on the vertices of a 
network reflecting the importance that each vertex plays in the network. Thus, for the 
case of weighted directed networks, one is looking for maps 

wdigraphs —?■ rankings on vertices 

sending a weighted directed graph to a ranking on its set of vertices. Asking for a lin¬ 
ear order on vertices is clearly too much since there could be vertices playing equally 
important roles in the network. A ranking on a set X is a pre-order dehned by a map 
r : X —)■ M such that x < y ii and only if r{x) < r{y). Clearly, different maps r may 
give rise to the same pre-order, i.e. to the same ranking. 

Finding a suitable map as above is not an easy task, and it is probably a problem 
with no universal solution. A hrst approach to this problem is via the total degree map 

G : wdigraphs —> rankings on vertices 

which orders vertices according to their total degree, i.e. the sum of weights of edges 
reaching or leaving a given vertex. With this ordering strongly connected vertices are 
deemed as having greater importance for the network. 

We work with networks of influences. As it often happens in mathematics, it is better 
to leave some terms undehned and let examples tell us the intended meaning. So, allow 
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US to review a few of examples of networks of influences. 


Author Citation Networks. 

Vertices in these networks are authors of scientihc publications. An author j has 
exerted an influence on author if there is at least one publication of i citing a publication 
of j. The matrix of direct influences is given by 

^ (j publications of i citing a publication of j 
(j publications of i 

The PWP method is designed to study indirect influences in networks of this sort. 

Business Providers Networks. 

Vertices are businesses in a group or economic sector. A business j influences business 
i if j provides products (or services) to i. The matrix of direct influences is given by 

^ amount spend by bussines i buying products from bussines j 

budget of bussines i 

International Trade Network. 

Vertices are countries. A country j exerts an influence on the economy of country i 
if there is trade between i and j. Thus the network itself is undirected, but the weight 
on edges do take direction into account. The matrix of direct influences is given by 

n.. = + Ejj 

Q ’ 

where Jjj is the total amount that i imports from j, Eij is the total amount that i 
exports to j, and Ci counts the total amount of international trade of country i (total 
of imports plus total of exports.) This example have been studied by Diaz and Gomez [6]. 

Process-Matter Networks. 

In these networks we have two types of vertices: black vertices for processes and 
white vertices for matter. The network represents a production system consisting of sev¬ 
eral processes. Each process takes some materials as input, and produces other materials 
as output, which may in turn be the inputs of other processes, and so on ... A matter 
vertex m influences a process vertex p if m is one of the inputs that p needs to operate; 
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a process vertex p influences a matter vertex m, if m is one of the outputs that p produces. 

Clearly these type of networks can be use to model a host of phenomena: chemical 
networks, metabolic networks, business organization systems, etc. The choice of weights 
for a process-matter network depend on the intended application, but they will always 
have the form: 

f 0 M \ 

[P 0 ) 

where indices for the hrst block of columns represent processes, indices in the second 
block of columns represent materials, the matrix P represents the influences of processes 
on materials, and the matrix N represents the influences of materials on processes. 

It is often useful to study influences among processes themselves, and also influences 
among materials themselves. This can be readily achieved by using, respectively, the 
product matrices 

MP and PM. 

In words, a process p influences a process q if an output of p is used by q as input, and 
similarly a material m influences a material n is there is at least one process that uses 
m as input and produces n as output. Figure [U shows an example of a process-matter 
network. Figure [2] displays the associated networks of influences among processes them¬ 
selves, and among materials themselves. 

With these examples in mind we go back to our main problem. From the viewpoint 
of networks of influences, we say that the total degree map G dehnes a ranking of vertices 
based on their direct importance. Although of interest, ranking by total direct degree 
fails to acknowledge that some direct influences may be short lived (if exerted over an 
isolated vertex) or may be enhanced (if exerted over a highly connected vertex.) Thus 
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Figure 2: Left: Associated Processes Network. Right: Associated Matter Network. 

although G should not be overlooked, it should not be regarded as the unique or final 
answer. 

A deeper approach to our problem should take indirect influences into account; how¬ 
ever there are many alternative ways for doing so. To proceed forward we look for a 
map 

R : wdigraphs —> rankings on vertices 

of the form: 


wdigraphs —> wdigraphs —> rankings on vertices 

where the map T transforms a network of direct influences into a network of indirect 
influences, and L is any ranking method computed from direct influences. 

To be able to pick a particular map T among the many possibilities some choices 
must be made. The PWP method for counting indirect influences [5] was founded over 
the following principles: 

1. Indirect influences arise from the concatenation of direct influences. 

2. Indirect influences do not arise in any other way. 

3. The weight of a concatenation of direct influences is proportional to the product 
of the weight of the direct influences that it comprehends, suitable modified to be 
compatible with our next principles. 

4. As a rule, the longer a concatenation of direct influences, the lesser the indirect 
influence exerted by it. 

5. Stochastic direct influences should generate stochastic indirect influences. 
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6. Robust convergency properties are expected. 

7. The map should be equivariant under to conjugation. 

The PWP method is dehned via the mapping 

T 

wdigraphs —> wdigraphs, 
which can be described in simpler terms as a map 

M^M) ^ M^M), 

from real square matrices of size n to itself, since we can identify (simple) weighted di¬ 
rected graphs with their adjacency matrices once a linear order has been hxed on their 
vertices. 


The map T depends on a real parameter A > 0, and is given on D G 

/ Er=i 


(R)by 






T(D,A) = 


— 1 


E oo 
k=l 


^ 

fc! 


The map T satishes our six requirements. Properties 1, 2, and 3 hold since the entries of 
matrix count weighted-paths of length k, i.e. concatenations of k direct influences. 
The dividing factor T implies property 4. Property 5 justihes the inclusion of the nor¬ 
malizing factor e^. Property 6 is insured by the factors Property 7 can be easily 
verihed. 


There is a number of other good choices for the map T, for example: the Katz index 
j9] , the MICMAC of Godet [8] , the PageRank of Google [H El [10] , the Heat Kernel of 
Ghung |1], and the communicability method of Estrada and Hatano [7]. These well-tested 
methods fail to satisfy some of the above requirements for the following reasons: 

• The Katz index (1953) when suitable extended for weighted networks and normal¬ 
ized, satishes most of our requirements; however its convergency properties are not 
as strong as desired, this is the main reason why the factors ^ are included in the 
PWP method. 

• MIGMAG (70’s) considers paths of a hxed length k, thus it fails to satisfy 1. 

• PageRank (1999) computes indirect inhuences by hrst transforming the matrix of 
direct inhuences into a Markovian matrix D; in this process inhuences are created 
using a mechanism diherent to concatenation. Thus PageRank fails to satisfy 2. 
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• Heat Kernel (2007) comes pretty close to satisfying our principles. However, it 
introduces indirect self-influences not coming from the matrix D, thus it fails to 
satisfy property 2. This is the main reason why the PWP method uses the function 
e^, instead of the exponential map e*. 

• Communicability (2007) corresponds (after normalization) to the A = 1 case in the 
heat kernel method; so the same arguments as above apply. Note that in both 
cases, starting with vanishing direct influences one obtains non-vanishing indirect 
influences. In contrast, in the absence of direct influences, the PWP method yields 
vanishing indirect influences, as also does the Katz index. 

Pre-composing any ranking based on the matrix of direct influences with the PWP 
map one obtains a new ranking that takes indirect influences into account. We consider 
three ranking methods: 

• The ranking by total degree G, which after composition with the PWP map we call 
ranking by importance. 

• The ranking by outgoing degree F, which after composition with the PWP map we 
call ranking by indirect influence. 

• The ranking by incoming degree E, which after composition with the PWP map 
we call ranking by indirect dependence. 

The orderings by incoming, outgoing, and total degrees may be thought, respectively, 
as maps 

E, F, G : M„(M) —)■ rankings on [n] 

where for a matrix D the place of a vertex in the respective orders is proportional to the 
values of the maps 

n n n 

Ei = Fi = and G* = ^(Aj + A*)- 

j=l j=l j=i 

Pre-composing these maps with the PWP map we obtain, respectively, maps 

E, F, I : Mji(M) — y rankings on [n] 

where for a matrix of direct influences D the place of a vertex in the respective orders is 
proportional to the values of the maps 

n n n 

Ei = ^A, Fi = ^A, and k = '^{Tij + Tji), 

j=l j=l j=i 
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defined in terms of the PWP matrix of indirect inflnences T = TD. 


Since there is only a hnite nnmber of rankings on a hnite set, the maps E, F and I 
can be continuous only if they are constant. Our main interest in this work is to study 
the regions of continuity of these maps, i.e. the regions where they are constant maps. 
Note that the rankings imposed by E, F and / on the set of vertices depend ultimately 
on D and A. We consider continuity with respect to D and A separately, as they play 
quite different roles: discontinuity with respect to D is taken as a manifestation of data 
sensitivity, in itself a positive phenomena, whereas discontinuity with respect to A is 
taken as a sign of the care required in the choice of A in the applications. The main 
question is whether or not changing A will let the PWP method to impose a more or 
less arbitrary ranking among the vertices. Our results, although partial, indicate that 
the opposite is the case: even if several rankings may result as A varies, they are seldom 
arbitrary, and evolve following a rigid pattern. 

2 Data Sensitivity of the PWP Method 

In this section we analyze the continuity of the map 

/ : M„(M) —)■ rankings on [n] 

obtained from the ordering by total degree after an application of the PWP map T{D, A). 
We set A = 1 and consider discontinuities of the map I with respect to matrix D of direct 
influences. Note that the matrix D is provided by the user, and we are going to test the 
sensitivity of the method to small changes in the matrix D. As expected, small changes 
in D can lead to quite different rankings on [n], stressing the need for a judicious choice 
of data. 

Let Lg be the linear graph with 6 nodes, see Figure 01 and let L6(e) be a perturbation 
of it obtained by adding a new edge of weight e as shown in Figure [HI For e = 0 we 
recover the graph Lg, further studied in Section [3], and as e grows the new edge becomes 
more relevant and modihes the ranking by importance of the vertices. 

Constructed numerically. Table [2] shows how the ranking of vertices by importance 
changes as e varies. Already for e = 0.01 we see that a change occurs, for example, the 
pairs of vertices {3,4}, {2, 5} and {1,6} have the same importance for e = 0, but for 
e = 0.01 all of them have different importance. Further changes in ranking were located 
at e = 0.28, 0.69, 2.1, 2.8, 7, and 23.9. Note that the vertex 4 is on the leading position 
up to e = 2.8, where the vertex 2 becomes the most important one. Note also that 
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Figure 3: Left: Graph L6(e). Right: Graph C^^e). 

vertex 3 begins at the leading position, and for e > 7 becomes the less important ver¬ 
tex. After we reach 23.9 our numerical experiments showed no further changes in ranking. 


This fairly simple example, already shows the high data sensitivity of the ranking of 
vertices by importance based on the PWP method. 
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Let us consider a second example, the circuit graph C^i^e) with 6 nodes extended by 
a new edge of weight e as shown in Figure [3l For e = 0 all vertices, as will be shown in 
Section m are equally important. But as soon as e reaches 0.0001 we hnd the ranking 


3,6 > 1,2 > 4,5, 

which, according to our numerical calculations, remains stable for higher values of e. This 
example displays both behaviours a highly sensitive one at the beginning, followed by a 
fairly stable one for higher values of e. 


We close this section with a remark on the impact that a change of scale, at data-level, 
have on the applications of the PWP method, i.e. we let our matrix of direct influences 
D be replaced by a new matrix of direct influences cD, with c G M+. We have that: 


T{cD,X) = 


1 °° \k 


k=l 


k\ 


1 

k\ 

+ k=l 


ocA 


-T{D,cX). 


k\ 





















As ^ is a positive real number, it does not affect the rankings by dependence, influence 
or importance. Therefore in the applications of the PWP method rescaling data by a 
factor amounts to rescaling A by the same factor, and thus stability with respect to data 
rescaling is a particular case of stability with respect to changes in A. 

3 Stability on Linear Graphs 

In this section we begin our study of the stability of the PWP method with respect to 
changes in the parameter A > 0. Recall that whereas the matrix of direct influences D 
comes directly from the user’s knowledge and experience, the parameter A comes from 
the PWP method itself. So, it is important to have a good control of the dependence of 
the PWP method on the choice of A. 

We consider the linear directed graph L„ with n vertices. Before considering the 
general case we deal with three simple but illustrative examples, namely, the graphs 
L 2 , L 3 , and Lg shown in Figure 01 





Figure 4: Linear graphs L 2 , L 3 and Lg. 


For L 2 , it is easy to check that for all values of A both vertices 1 and 2 have the same 
importance given by 


A 


— 1 


n=0 


A” 


n\ 


where Bn are the Bernoulli numbers |5]. Thus as A goes to inhnity the importance of 
both vertices goes to zero, meaning that the network becomes less connected, each vertex 
becomes more isolated, and thus its importance for the network decreases. See Figure [5l 


Regarding influences, vertex 1 has influence while vertex 2 has influence 0. Thus 
vertex 1 is always on top of vertex 2 in the order of influences, however as A goes to 
inhnity both inhuences tend to be equal to zero. Therefore for the graph L 2 the PWP 
method is fully stable both in importance and in inhuence. 
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Figure 5: Curves of importance for the graphs L 2 , L 3 , and Lg. 


Let us now consider the case of the linear graph L 3 with three vertices. Regarding 
importance, one can show that there is symmetry around the center of mass of the 
graph, as we will show below happens for all linear graphs L„. Thus vertices 1 and 3 
have the same importance, and therefore we only need to consider vertices 1 and 2. Here 
something interesting happens, as shown in Figure |5l for small values of A vertex 2 is the 
most important, but as A grows vertex 1 overcomes vertex 2 in importance. The matrix 
of indirect influences T for L 3 is given by: 


T = 


And thus the importance of vertex 1 and 2 are given, respectively, by: 



h{\) 


2A +A^ 

2e^ 


and 


/2(A) = 


2A 



Thus /i( 2 ) = 12 ( 2 ), Ji(A) < / 2 (A) for 0 < A < 2 , and h^X) < /i(A) for 2 < A. 
Therefore, for the graph L 3 the three possible rankings of its vertices by importance, 
after taking symmetry into account, actually occur, see Figured Influences in turn are 
fully stable as we have that 


m) 



R2(A). 


Let us consider the graph Lg 
given by 



with six vertices. 
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31 
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Ai 
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0 
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The matrix of indirect influences is 

0 0 \ 

0 0 
0 0 
0 0 
0 0 
A 0 y 
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By symmetry it is enough to consider the vertices 1, 2 and 3, with importance given by 


\2 \3 \4 \5 \2 \3 \4 \2 \3 

elh{\) = Y + '^ + ^ + ^’ e+/ 2 (A) = 2-^ +Y + -Y +Y’ ®+^3(A) = 2A + 2y+ —. 


.A 


Figure [5] shows the evolution of / 3 (A) as A varies. For A small we have the ranking 
1 < 2 < 3 in importance. Thus, initially vertex 1 is the less important one, and as A 
grows it hrst overcomes vertex 2 and then overcomes vertex 3 reaching the top position. 
Later on vertex 2 overcomes vertex 3, the ranking 3 < 2 < 1 is achieved, and it remains 
stable for large values of A. Indirect influences in turn are given by 


\2 \3 \4 \5 \2 \3 \4 \2 \3 

elFi(A) = A+y+ y+ y+ y, ^+^2 (A) = A + y + y + y , C+Fg ( A) = A + y + y , 




and thus the ranking 1 > 2 > 3 is stable for all values of A. Note however that the three 
values for importance approach 0 as A goes to inhnity. Thus, although the comparative 
values change in order, the overall values converge to zero. 


Next, we consider the linear graph with n vertices. The matrix D of direct influ¬ 
ences and the matrix T of indirect influences are given, respectively, by 


D^J = 


1 if i = j + 1 , 
0 otherwise, 


y-j 


if i > j. 


P^T- = 

J- ij 


0 otherwise. 

Proposition 1. The importance of the vertex j in the graph is given by 

j-i 


'-‘iW = E 


A^- 




X^-3 


1=1 ^ ^=]+l ^ ^ 

or equivalently: 

1. Forj = l,n, we have that e^/i(A) = e^J„(A) = A -|- ■ 

2. For 1 < i we have that 


2 A 2 


— 2A -h • • • -h 2 — 


A^' 


-1 


A^' 


(i -1)! ^ jF 


3. For ^ 


< j < n, we have that 


ellj{X) = 2A + ■■■ + 2 


Xn-j 


+ 


yn-j+l 


in-j)\ {n-j + l)\ 


+ 


A" 


(n-1)! • 


+ 


A«-i 
(n -i)!’ 


+ ■■■ + 


A^ 


/-I 
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4. The importance Ij{X) —)■ 0, as A —)■ cx). 


Proof. The hrst identity follows directly from the dehnitions. The other three identities 
follow from the hrst after specialization and simple changes of variables. Part 4 follows 
from the previous formulae. □ 

Next result gives us the symmetry in importance around the center of mass for the 
the graphs L„. 

Proposition 2. The importance of vertices in the linear graph is invariant under the 
change j —> n + 1 — j, i.e. we have for j G [n] that 


/,(A) = /„+i_^(A). 


Proof. We already know from Proposition [T] that /i(A) = Assume 1 < j < 

thus we have that n>n + l—j> and thus Proposition [T] implies that: 


--^n+l-j(A) — 2A + ■ • • + 


2 \n-{n-j+l) 


+ 


yi-{n-j+l)+l 


{n — {n — j + 1))! {n — {n — j + 1) + 1)\ 


+ ■■• + 


{n-j)\ 


2A^-i A^' 

^ ^ ^ 7 ! ^ 


\^~3 


□ 


By symmetry we only need to consider vertices to the left of the center of mass: 
namely for n = 2k or n = 2/c — 1 we only need to consider vertices j such that 1 < j < k. 


Consider the maps of importance Ij : (0, oo) —> M and the corresponding curves 

{ {X,Ij{X)) I A > 0 } C 

Our next result shows that the order in importance for large A is the reverse of the order 
in importance for small A. 

Lemma 3. For n = 2k or n = 2k — 1 we have that: 

1. If A is small enough, then Ii{X) < Ij{X) for 1 < i < j < k. 

2. If A is large enough, then /j(A) > Ij{X) for I < i < j < k. 
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Proof. According to Proposition [H for the smallest potency of A at which Ii{X) and Ij{X) 
differ are, respectively, of the form 

A* , A* 

— and 2—. 

^! ^! 

These terms control the behaviour of Ii{X) and Ij{X) for small A, and thus Ii{X) < Ij{X) 
since i < j. For large A, we look for the largest powers in A in A (A) and Ij{X) which are 
given, respectively, by 

An-i yn-i 

{n-z)\ (n-j)!' 

Since i < j, we have that /j(A) > Ij{X) for large A. □ 

Next we make a rather plausible statement which we have been able to verify numer¬ 
ically in many instances. 

Conjecture 4. Consider the linear graph L„, with n = 2k or n = 2k — 1, and let 
1 < * < J < The curves of importance /j(A) and Ij{X) intersect each other in a unique 
point A = Cij G (0, cx)). These intersection points occur in the order 

Cii < Cjm for i < I < k and j < m < k. 


The meaning of Conjecture 0] is that the curve A (A) begins at the bottom and crosses 
all other curves, hrst / 2 (A), then / 3 (A) and so on, until it reaches the top; next the curve 
/ 2 (A) raises from the bottom to the second highest position, just below A (A), crossing 
the curves / 3 (A), / 4 (A), etc, in exactly that order. After all crossings have taken place 
the reverse order to the original one has been achieved, i.e. /i(A) > / 2 (A) > ..., and this 
order remains stable up to inhnity. 

The numerical evidence for Conjecture 0] is quite solid. Figure [ 6 ] suggests that it holds 
for the linear graph Ln . 


Next result allow us to locate the crossing points of consecutive curves of importance. 


Theorem 5. Let L„ be the linear graph with n = 2k or n = 2k — l, and 1 < i < k — 1. 
The curves of importance /j(A) and /i+i(A) intersect each other in a unique point 
A = Ci^i+i e (0, 00 ) given by 


Q,2+i — 


{n — i)\ 


i\ 


1 

n — 2i 


Thus for n = 2k and n = 2k — 1 we, respectively, have that 


-' 1 , 1+1 


{2k - i)!\ 2 (^ 


i\ 


and 




(2k — i — 1 )!\ 2(k-i)-i 


ii 


13 










Figure 6: Curves of Importance for Ln. 


Proof. The crossing point is defined by the equation 


2A + + 

2A + + 


+ 2 ——— + — + 
(z —1)! ^! 


A* 

+ 2 - + 
z! 


yi+l 

(z + 1)! 


+ 


+ 


+ 


A" 


which after cancelling terms is equivalent to the equation 


An-i 


(n — 1 }\ 


A* 


Z! 


with a unique solution given by 


Q,i+1 


[n — Z !\ 


Z! 


(n — z)! 

An-i-l 

(n — 


1 )!’ 


□ 


Corollary 6. Let L„ be the linear graph with n = 2k or n = 2k—1, and let 1 < z < k—1. 
1. For n = 2k, k > 2, the crossing points ci ,2 and Ck-i,k are given by 
Cl ,2 = {2k-l)\^ and Ck-i,k = {k{k + 1))P 


2. For n = 2k — 1, k > 2, the crossing points Ci ,2 and Ck-i^k are given by 

Cl ,2 = {2k-2)\^i^ and Cfc-i,fc = k. 

The following result may be regarded as further evidence in favor of Conjecture HI 

Theorem 7. The intersection points Cj,j+i of the curves /i(A) and Ji+i(A) occur in the 
order 

Cl,2 < C2,3 < . < Ck-2,k-l < Ck-l,k- 
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Proof. We use a simple fact for positive integers: if c < a — 1, then 

ac < (a — l)(c + 1 ). 

Choose i such that 2{i + 1 ) < n, i.e. such that n — 2z — 2 > 0 or equivalently 

i + 1 < n — i — 1. 

We have the following chain of equivalent inequalities 


^ C^+1,^+2; 

t\ ) ^ ( (* + !)! 


(n-*)"-2'-2(z + l)-2*-2(, + i)!2 < 

(n - + 1 )’^-"'-' < (n-^-l) 2 •••(^ + 2 )^ 


((n-i)(z + l)) 


n—22 —2 


< {n — i — lY ■ ■ ■ {i + 2y 


Using the fact mentioned at the beggining , it is clear that in order to show the latter 
inequality it is enough to check that 


{n — i){i + l) < {n — i — l){i + 2), 

which holds since z + l<n — z — 1 . □ 

Example 8 . Set n = 12 and z = 3. In this case the inequality 

((^_,)(z + l))-2-2 < ^ri-i-lf---{i + 2f 
simply says that (9.4)"^ < (8.7. 6 .5)^, or equivalently 

(9.4)(9.4)(9.4)(9.4) < (8.5)(7.6)(7.6)(8.5). 

Note that if Conjecture 0] holds, then there are exactly 

k{k-l) 

2 

crossing points among the curves of importance /j(A), and thus that same number of 
different orderings by importance on the vertices of L„ as A varies, where n = 2k or 
n = 2k — 1. Indeed vertex 1 begins as the less important and have to surpass k — 1 
vertices to reach the top. Then vertex 2 have to surpass k — 2 vertices to reach the 
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second position, etc. Thus the number of crossing is the sum of the hrst k — 1 natural 
numbers, yielding the desired result. 


Although the number of reachable orderings grows to inhnity, it is nevertheless a 
negligible quantity, for large n, relative to the number of all possible orderings as 


him 

k^oo 


k{k-l) 

m 


0 . 


Therefore even though it is possible to reach many different orderings by choosing an 
appropriated A, a random ordering will not be reachable. 


According to Conjecture IU for the linear graph L„ with n = 2k or n = 2k — 1, the 
hrst change in the ordering of vertices by importance occurs at Ci^ 2 , while the last change 
occurs at Ck-i,k- The interval of stability (0, 01 , 2 ) is dominated by the direct inhuences, 
while the interval of stability {ck-i^k, 00 ) is dominated by long indirect inhuences, yielding 
the reverse ordering. All reordering happens in the interval [ci^ 2 -i Ck-i^k\- 


Finally, we consider the ranking of the vertices of L„ by inhuence, and in this case we 
hnd full stability. 


Proposition 9. For any A > 0, the ordering by indirect inhuences on the vertices of the 
linear graph L„ is given by 1 > 2 > • • • > n. 

Proof. It follows since Fn = 0 and for 1 < z < n we have that: 


Fi{X) — A + — + 


A*-^ A* 

+ j-. -^ + -r + 

(z — 1)! z! 


+ 


A" 


(rz — z)! 


□ 


Note however that while the ordering is completely stable, all inhuences approach 0 
as A goes to inhnity. 


4 Stability on Circuits 

In the previous section we saw that linear chains of direct inhuences tend to have a desta¬ 
bilizing ehect on the applications of the PWP method with respect to changes in A, to 
the point of allowing complete reversal in ordering for small and large values of A. In this 
section we are going to see that, in contrast, the presence of circuits have a stabilizing 
ehect with respect to changes in A. 
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Consider the circuit on the Z„ group of integers module n given by 


n 


0 . 


For i,j, G Z„ and k E [n], the matrices of direct and indirect influences are given by 


Dij 


1 if z = i + 1, 
0 otherwise, 


and 


( 2 ^ 7 ^ 


^k+ln 

^ {k + ln)V 


Note that Tj+fcj(A) does not depend on j, thus we can use the simpler notation Tk{\). 
Next result shows the full stability of the PWP for circuits: stability in the rankings by 
importance, in influence, and even in the relative strength of indirect influences for small 
values of A. 


Theorem 10. For i,j G Z„ we have that: 

1 . All vertices in Z„ have equal importance, that is /i(A) = /j(A). 

2. All vertices in Z„ have equal influences, that is -Fj(A) = Fj{X). 

3. For A G (0,2) indirect influences are ordered as follows: 

Ti(A) > ••• > T^A). 


Proof. Properties 1 and 2 follow from the identities 

F,(A) = Ti(A) + ■■■ + T,(A) and /,(A) = 2F,(A). 
Property 3 is shown as follows. Recall that for k E [n — 1] we have 


nw = E 


1=0 


{k + ln)\ 


thus the inequality T’fc(A) > Tfc_|_i(A) holds for each summand in the respective series 
expansions if and only if 

yk+ln yk+l+ln 


> 


or equivalently 

Thus for A < 2, we have that 


{k + ln)\ (A: + 1 + ln )\’ 
A < k + 1 + In. 


A < 2 = 1 + 1 < k + 1 < k + 1 + In, 

and we obtain the desired inequality. □ 

Thus we see that the presence of circuits in a complex network have a stabilizing 
effect in the applications of the PWP method with respect to changes in A. 


17 



5 Stability on M-Diagonalizable Networks 


Assume that our matrix of direct influences D is diagonalizable in R, i.e. there is an 
invertible matrix A G M„(R) such that D = AEA~^ where £" is a diagonal matrix whose 
entries Ea give the eigenvalues of D. This condition holds, for example, if the eigen¬ 
values of A are all real and distinct, a generic condition among matrices with only real 
eigenvalues. 

The PWP map is equivariant with respect to conjugation, thus we have that 

T(D,A) = AT{E,X)A-\ 

The matrix T(i?, A) is rather simple to compute, indeed it is a diagonal matrix with 
entries 

T(i?,A),, = 

e+ 

Theorem 11. For a R-diagonalizable network there can be only a hnite number of 
changes in the ranking of the vertices of the network by importance or by influence. 

Proof. We must show that only a hnite number of crossing points may occur, both for the 
curves of importance and the curves of indirect inhuences. Since D is a diagonalizable 
matrix in R it has n real eigenvalues (counted with multiplicity). Assume it has m dif¬ 
ferent eigenvalues which we write in increasing order di < ... < dm- The key observation 
is that the entries of the matrix 

e^T(D,A) = eXAP{E,X)A-^ 

are all of the form 

n J_ _L /-» 

+ * * * + 7 

and therefore the functions of importance e^Ij{X), and the functions of inhuence 
e^Ej{X) are also of the same form. Therefore, hnding the crossing points for the curves 
of importance or inhuence boils down to hnding solutions to equations of the form 

+ ■ ■ ■ + = 0 . 

Some of the resulting equations may be trivial, i.e. all the coefficients Oi may be zero, 
meaning that some of the importance or inhuences functions are identically equal, which 
actually reduces the scope of possibilities for crossing points. We are going to show that 
as soon as one of these equations is non-trivial there can only be a hnite number of 
solutions. Since — 1, the equation above is equivalent to a equation of the form 

aie‘^1'^ -h ■ ■ ■ + = a (E) 
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with a = ai + ... + am- Without lost of generality we assume that 7^ 0. If it is the 
only non-vanishing coefficient, then ([E]) reduces to = am which is either trivial if 

dm = 0, or has no solution at all. 

If another coefficient besides am is non-zero, we may assume without lost of generality 
that Oi 7^ 0. We show that if ([E]) has inhnitely many solutions they must be contained in 
a bounded interval around 0. Assume that dm > 0, otherwise all the eigenvalues di must 
be negative, and then as t —)■ 0 we have that 

aie‘^^^ ■ ■ ■ - 1 - ame’^^^ 0 , 

and therefore ([E]) has no solutions for large AifoT^O. Ifa = 0, then ([E]) is equivalent to 

-F ••• -F = -ai, 

and we are back in the case where the last coefficient dm — di is positive. 

Dividing ([E]) by and letting A go to inhnity we hnd that = 0, a contra¬ 
diction, and thus no large A can be a solution of ([E]). So, if ([E]) have inhnitely many 
solutions they must be contained in an interval around 0, and these solutions must have 
an accumulation point. Let us show that this point of accumulation can not be 0. We 
argue by contradiction. Assume that ([E]) has inhnitely many solutions ti with —)■ 0 as 
Z —)■ cxD. Then the derivative of the left-hand side of ([E]) will have inhnitely many zeroes 
accumulating at 0, and thus the second derivative will also have inhnitely many zeroes 
accumulating at 0, and so on ... Therefore the coefficients ai,... ,am must be such that 
the identities 

aidi ■ -1- amdm = 0, 

hold for Zc G N>o. Dividing the equation above by d^ and letting k go to inhnity we hnd 
that Om = 0, a contradiction. 

Finally, assume that our inhnitely many solutions ti of ([E|) have an accumulation 
point c 7^ 0. Set ti = c + si, then the points si give inhnitely many solutions to the 
equation 

accumulating at 0. Therefore we must have that 

= • • • = = 0, 


and so Oi = • • • = = 0. 


□ 
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Remark 12. If we allow complex eigenvalues, then besides exponential functions, trigono¬ 
metric functions may appear in the calculation of the matrix of indirect influences. 
Trigonometric functions can have inhnitely many crossing points. 

Thus for a diagonalizable network only a hnite number of changes in ordering either 
by importance or influence can occur, and thus for such networks there exists G (0, cxd) 
where the hrst change in order occurs, and a A* G (0, oo) where the last change in order 
occurs. Thus the PWP method for A in the interval (0, A^) is dominated by the direct 
influences, and the PWP method for A in the interval (Aj, oo) is dominated by the indirect 
influences. All reordering happens in the interval [A^, A*]. 

6 Conclusion 

In this work we have considered the stability of the PWP method for ranking vertices in 
a complex network by importance and indirect influences. We have found that the PWP 
method is quite sensitive with respect to data, a fact that we regard as being positive. 
Stability with respect to the parameter A seems to involve, at least, two opposite forces. 
On the one hand, long directed path increases instability, to the extreme of allowing full 
order reversal, and other hand the presence of circuits tends to stabilize the applications 
of the method, again to the extreme of allowing full uniformization both in importance 
and influence. Must networks, of course, include both directed paths and circuits, and the 
stability of the PWP method with respect to changes in A will involve a subtle balance 
between these opposite forces. We have shown that for a network diagonalizable in M only 
a hnite number of changes in ranking may occur, both in importance and in inhuence. In 
contrast, the presence of complex eigenvalues opens the door for inhnitely many changes 
in order. The next challenge, left for future research, is to study the stability of the PWP 
method for randomly generated networks. 
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